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Potential steps naturally develop in graphene near metallic contacts. We investigate the influence 
of these steps on the transport in graphene Field Effect Transistors. We give simple expressions to 
estimate the voltage-dependent contribution of the contacts to the total resistance and noise in the 
diffusive and ballistic regimes. 



Graphene's distinctive band structure gives rise to ex- 
citing new transport properties and promising applica- 
tions for carbon-based electronicsii^i^. When measur- 
ing the conductance or current noise in a nanotube or 
a sheet of graphene, the properties of the contacts can 
matter as much as the electronic structure of the nan- 
otube or graphene itself. In semiconducting nanotubes 
or graphene nanoribbons, it is known that Schottky bar- 
riers develop at the metallic contacts^. Charge trans- 
fer between a metal and a wide graphene strip induces 
potential steps whose shape may differ strongly from 
usual Schottky barriers due to the semimetallic and two- 
dimensional nature of graphene. The existence of such 
metal-induced potential steps was inferred experimen- 
tally from the transport properties of a graphene strip 
with various contact geometries^. More direct evidence 
for these steps comes from optical mapping of the poten- 
tial landscape across a graphene device^. Recent theo- 
retical work on graphene-metal interfaces has been per- 
formed within the atomistic tight-binding theory^i^. 

In this paper, the conductance and the shot noise of 
graphene Field Effect Transistors (gPETs) with extended 
contacts are derived using the Dirac Hamiltonian for 
graphene. Near a single contact, we assume that the 
Fermi energy in graphene varies monotonically over a 
characteristic length d, and we solve the corresponding 
scattering problem exactly. If the transport is ballistic 
between both contacts, we predict oscillations of the noise 
and conductance as the charge density is increased in the 
sheet. When the density exceeds the one under the con- 
tact, the noise minima might be zero and correspond to 
perfect transmission between the contacts. Such realiza- 
tions of a noiseless gFET are caused by Fabry-Perot res- 
onances and require low doping by the contacts. In the 
diffusive regime, we show how the total resistance and 
Fano Factor of the whole gFET depend upon the contact 
resistance and Fano factor of each contact, which is rel- 
evant in interpreting recent experiments on shot noise in 
non-suspended graphenei2i>ii. 

Before analyzing the gFETs' properties, it is useful 
to investigate transport across a single graphene-metal 
contact. We thus consider that a metal electrode cov- 
ers the left half-plane {x < 0) of an infinite graphene 
layer. We assume that the metal coating simply shifts 



metal/graphene interface as a function of transverse momen- 
tum ky for several lengths d of the potential rise. The ra- 
tio between Fermi wavevectors in bulk graphene and below 
the metal is set to fc^'/fc^^ — 3 (unipolar contact). A 
schematic of the device and corresponding Fermi wavevec- 
tor profile kp{x) are shown in the inset, b) same curves for a 

(ff) _ q 



bipolar contact with kp^'^ /k^j^ 



the Fermi level of the graphene underneath while pre- 
serving its pristine Dirac conesi^. Far from the contact, 
the type (n or p) and density of charge carriers in the 
right half-plane {x > 0) are tuned by a distant metal- 
lic gate. A continuous Fermi wave vector profile kpix) 
must therefore develop near the contact edge to match 
the asymptotic values kp{— 00) — fc^-* under the metal. 



and fci?(+oo) 



_ iJg) 



in bare graphene. The dynamics 



of the massless fermions can be safely described by the 
single- valley two-dimensional Dirac Hamiltonian, 



H = —ihvp [C'xdx + <7ydy — ikp{x)l\ 



(1) 



which is valid if the potential step is smooth on the scale 
of the lattice constant. Here vp is the Fermi velocity and 



kpix) = 4™^ + 



k 



(9) 



(2) 



our simple model for the space-dependent Fermi wave 
vector (Fig. [1] inset). The Pauli matrices {<Jx,cFy,az) op- 
erate on spinors ^'(a;,?/) whose components are the elec- 
tron amplitudes associated with each sublattice of the 
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honeycomb carbon crystal. It is worthy to note that a 
complete treatment of nonlinear screening and disorder 
near the contact should lead to a more complicated pro- 
file in the bipolar n-p cas&i^ii^. Nevertheless this can 
be accounted for by noting that our phenomenological d 



will in fact depend on /c^*' and kp' in a way that may 
ultimately be calculated. 

We now proceed to the derivation of simple formulas 
relating the transport properties (as functions of kp"^) 

to the parameters d and fc^'* characterizing the con- 
tacts. Here we consider a wide enough graphene strip 
to neglect edge effects^. Then translational invariance 
parallel to the junction (along the y axis) implies that 
the transverse wavevector ky is conserved, and all spinor 
wavefunctions take the form ^{x, y) = $(a;)e*'^''^ where 
<I>(x) = "^($i(x), $2(2;)). Aiming to determine the low 
temperature transport, we solve the Dirac equation at 
the Fermi level {H'^{x, y) — 0) which reduces to the one- 
dimensional equation 



dx^{x) = {kyaz +ikFix)ax)^{x). 



(3) 



This equation can be decoupled using the symmet- 
ric/antisymmetric combinations of components f±{x) — 
^i{x)zt^2{x) which obey the following scalar differential 
equations: 



iakp{x))fa = 0, 



(4) 



where a = ±. In the asymptotic regions |a;| 3> d, the 
solutions are plane waves fa{x) = e***''==^ which can be 
either exponentially damped or oscillatory depending on 
the sign of fc^ — kp — ky. In order to find the transmission 
across the potential step Eq. we now construct a 
scattering state containing a single oscillatory outgoing 
charge state (with \ky\ < {kp^'D in the region x — > +00, 
namely 



Here /ci™^ — s„ 



fc^ and Sm = sign(fc^'') indi- 



cates whether graphene is n- or p-doped underneath the 
metal. The Euler Gamma function is denoted T{z) and 
we have introduced the momenta k^"^ = kp^ — fc^'"'' + 

pfci"'' + <yk^\ with p,a ^ ±1. The corresponding reflec- 
tion probability is simply given by 



91: 



step 



r(i - iK++d)r{iK—d) 



r{l -iK+-d)T{iK-+d) 



(9) 



when all the waves are propagating, namely for \ky\ < 

min(|fc^'|, Finally, a remarkably simple formula 

is obtained for the reflection coefflcient of Dirac fermions 
across the potential step Eq. ([2]): 



step 



sinh(7rc?K+ ) sinh(7rc?K 
sinh(7rc?K++) smh^Trdu' 



(10) 



This expression is valid for any d and step po- 
larity, matching the known limits for transport 
across smootbi^i and abrupt {d 0) steps'^^, 

and interpolating between those limits. The ex- 
pression ([Toll is reminiscent of the reflection coef- 
ficient of a non-relativistic massive scalar particle. 



d\ — sinh ( nd 



(9) 



I sinh ^TTc 



1.(9) 



obtained by solving the Schrodinger equation in a similar 
potential landscape^"^. The richer structure of Eq. (fTU|) is 
associated with the Dirac nature of carriers in graphene. 
In particular, it indicates the absence of backscatter- 
ing at normal incidence (fci™'^'' = /c^'^^) for any height 
and width of the potential step, which is related to the 
orthogonality of incoming and reflected spinor states^^. 
As with standard impedance matching, the transmission 



e - at a; 



(5) 



Here k^^ 



{kfy 



k"^ is the longitudinal mo- 



mentum whose direction depends on the band index 
Sg = sign(fc^-') far on the right side. In Appendix A, 
we show that this asymptotic condition completely de- 
termines the solution of the Dirac equation H'i>{x, y) ^ 
on the whole x-axis. In particular, on the left side the 
wave consists of a superposition 



Ux) - fi'"'Kx) + frf\x) at X ^ 
of an incoming 

^ r(i-2zfci^^d)r(-2»fci"')d) (^,, 

^' T{iK—d)T{l-in++d) 
and a reflected charge carrier 

r(i - 2ik'iU)T{2ik'i^U) 



J cx 



{ref) 



V{iK-+d)T{l 



d) 



(6) 



(7) 



(8) 
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FIG. 2: Contact conductance (left) and Fano factor (right) as- 
sociated with a single potential step as functions of the Fermi 
wavevector in bulk graphene fc^' for several values of the 



(9)/l.('») 



dimensionless parameter fc^'d. 



The maximal value of the 



conductance is given by Go = 4e N/h where A'^ = kp W/n 



T^tcp = 1 — 9^stop of the unipolar contact tends toward 
unity when the distance d is increased (Fig. [1]) . In con- 
trast, for bipolar steps, the fcy-dependent transmission 
T^top goes from a broad curve at small c? to a sharp peak 
around fcj, = at large d (Fig. [T|) . Besides, the bipolar 
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transmission counter-intuitively increases when the po- 
tential barrier height is increased. Although the poten- 
tial step originates from charge transfer between a metal 
and graphene, these scattering properties are similar to 
the relativistic Klein tunneling^*' for which evidence is 
mounting in the context of transport through potential 
barriers created by local gate o-'^^i^-'^'^^i^^'^'^i^^'^^ . 

From the transmission probability T^tepj the conduc- 
tance and the Fano factor of a single contact are given 
respectively by 



4e2 W 
h '2tt 



(11) 



and 



F CA-'^h 



dfcy7^tep(l ^top)/ 



where W is the width of the graphene strip along y 



and k 



is related to the asymptotic density n'^' 
.(ff) 



at 



a; ^ d by the relation k'^p' = vthtJ^. These quanti- 
ties are strongly sensitive to the nature and density of 
charge carriers (Fig. The bipolar contact is clearly 
more resistive and noisier than the unipolar one, and this 
unipolar/bipolar asymmetry becomes more pronounced 
for smoother contacts (Fig. In the limit of very 



smooth potential steps {k'^^^'d — > oo), the Fano factor 
vanishes in the unipolar case whereas it saturates to a 
finite value F^onikp^ < 0) = 1 - 2^1/2 ^-j^g bipolar 
case (in agreement witbii). As expected a single con- 
tact becomes noiseless when the potential step vanishes 
(fc^^ — k'^p'^). In the sequel of the paper, we use the 
result Eq. PH)) to investigate the effect of the contact po- 
tential steps^'^ on the conductance and noise properties 
of gFETs. 

Recently, suspended gFETs have been achieved, re- 
sulting in increased mobility of the two-dimensional elec- 
tron gas22iS^i2^. In such devices electronic motion might 
become ballistic between source and drain. The whole 
structure can be described by two symmetric steps sim- 
ilar to Eq. and separated by a distance L, the so- 
called Wood-Saxon potential Eq. (jBl[) . The single chan- 
nel transmission exhibits Fabry-Perot-like resonances as 
the gate voltage varies. We calculate the transmission 
probability using a mapping between the problem of 
massless Dirac fermions in graphene and the one of mas- 
sive Dirac fermions in the one-dimensional Wood-Saxon 
potentiai^. The contrast of the whole interference pat- 
tern is controlled by the reflection probability 9'tstop(fcy) 
which depends strongly on the distance d, as discussed 
above, see Eq. PH)) . 
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In the large gate voltage regime, 

wavefunctions are oscillatory for all transverse channels 
and the corresponding transmission across the whole de- 
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FIG. 3: Conductance G of a ballistic Fabry-Perot cavity as 
a function of the dimensionless parameter k^p^L. The maxi- 
mal conductance (Go = Ae^N/h) is controlled by the number 
of available propagating modes in the leads A'' = fc^'iy/7r. 



Here k 



1.5, which means L = 27 nm for a density 
and the aspect ratio W/L is large (in 



practice W/L > 4 is sufficient ). For definiteness, we have 
assumed that the graphene is n-doped underneath the metal 
(fc^™' > 0). In this case, the minimal conductivity is reached 
at a nonzero negative density in the central channel. For other 
metals that produce p-doped graphene after charge transfer!, 
one should obtain the symmetrical curves (with respect to ver- 
tical axis k^^^L = 0) of the ones presented here, and the shift 
of the minimal conductivity would occur at positive density. 



vice can be approximated by 



^Stcpiky) 



%Mky) ~ 1 + 



(1 - mstcpiky))' 



■ sm 



iki'^L)) (12) 



when neglecting the corrections of order d to the effective 
size of the cavity. The shape of the cavity enters this for- 
mula through the reflection probability y^stcp{ky) ■ The 
actual position of the conductance peaks (Fig. [3]) is mod- 
ified because the effective width of the cavity changes 
with d. In this regime, the conductance of realistic cavi- 
ties is globally smaller (Fig. [3]) than the prediction of the 
square well-model (d = 0) with infinite doping below the 
electrodes (fc^'' = oo)^. Indeed the maximal conduc- 
tance {Ae^N/h) is controlled by the number of available 
propagating modes in the leads N — k^^^W/i: as soon 
as the doping of the two-dimensional electron gas in the 
middle part of the gFET exceeds that underneath the 
source and drain contacts. At large absolute gate volt- 



ages {kp^ < —k]P'), the bipolar conductance exhibits 
strong oscillations whose contrast increases as the steps 
become smoother (Fig. |3]). In contrast at kp'' > 
the unipolar cavity becomes fully transparent in the large 
d limit wherein the oscillations are lost. Similar fringes 
have already been observed recently in short Fabry-Perot 
devices created by local gating^S. We suggest making 
contacts with a metal which dopes the graphene very 
lightly in order to observe the Fabry-Perot interferences 
with the highest resolution in bipolar cavities. 
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We now discuss the more usual low gate voltage regime 



da) 



< 



wherein the conductance is determined 

by both evanescent and propagating modes. Besides be- 
ing nonuniversal, the minimal conductance is reached at 
nonzero carrier density ^ 0. Up to now, the mini- 
mum conductance was predicted to occur at a non-zero 
gate voltage (due to charge impurities) but at zero av- 
erage density^. Here, we predict that the presence of 
metallic contacts can shift this minimum to non-zero 
charge density, even without impurities. Besides, the 
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FIG. 4: Fano factor for a ballistic graphene Fabry-Perot cav- 
ity as a function of the dimensionless parameter kp^'L for 
various values of d/L. Same parameters as in Fig. 3, in par- 
ticular k^™^ L = 1.5 while the aspect ratio W/L is large^, 
namely W/L > 4. 

Fabry-Perot interferences should also induce signatures 
on the noise of clean gPETs. Here we demonstrate that 
the Fano factor of a ballistic gFET, 



(13) 



should exhibit spectacular oscillations as a function of the 
dimensionless parameter kp^L (Fig.|4]) for cavities oper- 



ated at gate voltages yielding 



for kp ^ > kp'^' , the Fano factor shows successive nodes at 
the gate voltages corresponding to peaks of optimal con- 



(9) 



> /c^*' . In particular 



ductance (Fig. [3]). On the bipolar side kp' < —kp', the 
nodes are replaced by local minima of -Fbaii corresponding 
to local maxima of the conductance. These oscillations 
have not been investigated either experimentally or theo- 
retically. Early theoretical works assumed infinite doping 
by the contacts^^,. 

We now consider the situation, relevant to non- 
suspended gFETs, where the carrier motion is diffusive 
between the source and the drain. Assuming phase- 
incoherent transport, the total resistance of the gFET 

is simply given by i?dif = Rshcct{k';^^ ) + 2Rcon{k';^^) and 
the Fano factor by 



36 



Fdif(4^^) 



••sheet ^shcct 



(2i?c 



-Rshcct)^ 



(14) 



where i?shcct and Fgheet are respectively the resistance 
and the Fano factor of the sheet, and i?con = Gcon 
the single contact resistance. Theory of the anomalous 
diffusion through the electron-hole puddles sea formed 
in graphene predicts a universal scale-independent Fano 
factor Fshcct — 1/3"^^ in agreement with a recent 
experiment^". In addition numerical studies also indicate 
a nearly density independent Fano factor (not universally 
equal to 1/3) from moderate to strong disorder— i^i. Ac- 
cording to Eq. (fT4|) . one expects -Fdif(fc^^) = Fgh, 



1/3 only when i?shoGt > Rcon 
in all devices measured in 



which might be the case 
Ref.— . In contrast, when 
-Rshcot < -Rcon, Fdifikp') depends more on i^con which 

. Interestingly 
decrease from 



makes it decrease as a function of 
another experiment reports on such a 



Fdi{{kp^ =0) — 1/3 to a lower value at large densitiesii. 
Nevertheless it should be emphasized that the incoher- 
ent theory underlying Eq. p^ is not valid for the very 
short graphene devices investigated in Ref.^-'^. It is thus 
necessary to consider the nonlocal transport properties 
of the whole device. Unfortunately the precise profile of 
the potential in this experiment is not known and prob- 
ably corresponds to L ~ c? making neither our study 
of the ballistic Fabry-Perot (Fig. [4]) nor the square well 
model of Ref.>32. quantitatively relevant. Nevertheless the 
qualitative behavior of these experiments is captured by 
these models: intrinsic noise close to the Dirac point and 
contact-dominated noise at larger densities. 

In conclusion, we have considered Dirac fermions scat- 
tering from one or two potential steps having each a char- 
acteristic length d. Such steps introduce additional dis- 
sipation localized at the source and drain of gFETs and 
also modify drastically the noise properties of such de- 
vices. In the ballistic regime, we predict that the presence 
of metallic contacts can shift the conductance minimum 
to a non-zero charge density, which is negative (resp. pos- 
itive) for a metal which dopes the graphene with elec- 
trons (resp. with holes). In addition we also suggest 
performing conductance and noise measurements on a 
suspended bipolar graphene Fabry-Perot structure22i22ii2S 
with low doping at the electrodes in order to observe en- 
hanced Fabry-Perot oscillations. Finally, the asymmetry 
between electron and hole transport is enhanced when 
these potential steps rise over longer distances. 



APPENDIX A: TRANSMISSION ACROSS A 
POTENTIAL STEP 



We consider the differential equation Eq.(l4|): 
rj, + ikUx)-kl-aikpix))U = 0, 



where 



(Al) 



(A2) 
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We introduce the new independent variable ^ = — e~^/'^ 
and seek for the solutions in the general form fa{x) — 
'^^ct(C)- The function WaiC) satisfies the Gauss 
hypergeometrical equation: 

e(l - OK + [c-{a + h+ < - ahw^ = 0, (A3) 

if the exponents are chosen as 



= -iki3)d and \ = ^i{k^f -¥p)d. (A4) 



Introducing v — ik!""^^ d, the parameters a, 6, c are given 
by 



a+h+l = 
a = 
h = 



1 + 2/i, 
1 + 2/i - 2A, 
- A + v), 
(pi - A - j^). 



(A5) 
(A6) 
(A7) 
(A8) 



In the region x +oo, namely ^ — > 0, the two 
independent solutions of the hypergeometric equation 
are Wa{0 — F{a,b,c;£_) and Wa{£,) — ^^~^F(a — c + 
1,5 — c + 1,2 — c;^)"^^. The corresponding functions 

fa{x) — ^'^WaiO the plane waves ^'^ — e*'^^"''' and 
= e-«'=i'''<i. We now construct a scattering state con- 
taining a single oscillatory outgoing wave in the region 
X +00, namely 



fa{x) 



at X ^ +00, 



(A9) 



where fci^"* 



''sV i^F^)^ ^ the longitudinal mo- 

mentum whose direction depends on the band index 
Sg — sign(fc^-') far on the right side. Note that although 

the sign of fci^"* can be either positive or negative de- 
pending on doping, the group velocity always describes 
a right-moving charge. 

From the general relational between the hypergeomet- 
ric functions of respective arguments ^ and 1/^, one can 
extract the asympotic behavior of w(^) in the region for 

^ — > —GO 



WaiO 



T{c)T{b~a) 
r(6)r(c- a) 

r(c)r(a- 



+ 



T{a)T{c-b) 



(AlO) 
(All) 



Consequently the structure of the wave fa{x) : 
5)~^w(^) in the region x 

(namely ^ — * — oo) consists in two parts 



-oo 



Ux) ^ fi'-'Hx) + fi^^f\x) at 
which are respectively the incident wave 



— oo. 



(A12) 



(A13) 



and the reflected wave 



r(a)r(c- 6) 



We have used 



-X I d \ —V 



-6+/J-A _ 



x/d\iy 



-ik 



(™), 



(A14) 

(A15) 
(A16) 



We emphasize again that e*'^^'"'^ is always the right- 
moving incident wave, although fci™'' can be either pos- 
itive or negative (because the projection of the group 
velocity is positive in both n-type and p-type doped 
graphene). Therefore the reflection probability is: 



9^, 



■step 



r(6)r(c- a) 



r(a)r(c- 6) 
r(^ - A - j/)r(i + A + /i - i^) 



(A17) 



(A18) 



r(/i - A + j/)r(i + A + /i + jy) 

which leads to Eqs. (|9ll0p in the text. 



APPENDIX B: WOOD-SAXON POTENTIAL 

The Wood-Saxon potential corresponds to two sym- 
metric steps 



kp{x) 



"F ' 

i-x) 



4'"' 1 X 



e{x) 



p-{x+L/2)/d ^ I p(x-L/2)/d I 



(Bl) 



We have checked that the transmission probability is 
given by the formula of Refj^, namely 7baii(fcy) — 



A 



C 



(m + - A- 

r2(-2/x)r2(-2i/) 



A) 



A2 



A2 



A) 



A)- 



(B2) 
(B3) 
(B4) 



We denote S(a, b) = T{a)T{b)/T{a + b) the beta function. 



and /I 
in appendix A 



-i{k 



)d like 
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